Local Measures of Entanglement and Critical Exponents 
at Quantum Phase Transitions 
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We discuss on general grounds some local indicators of entanglement, that have been proposed recently for the 
study and classification of quantum phase transitions. In particular, we focus on the capability of entanglement 
in detecting quantum critical points and related exponents. We show that the singularities observed in all local 
measures of entanglement are a consequence of the scaling hypothesis. In particular, as every non- trivial local 
observable is expected to be singular at criticality, we single out the most relevant one (in the renormalization 
group sense) as the best-suited for finite-size scaling analysis. The proposed method is checked on a couple 
of one-dimensional spin systems. The present analysis shows that the singular behaviour of local measures of 
entanglement is fully encompassed in the usual statistical mechanics framework. 

PACS numbers: 05.70Jk, 03.75.Ud, 75.10.Pq 



The entanglement properties of condensed matter systems 
have been recently object of intensive studies 1 1, 2], especially 
close to quantum critical points (QCP) where quantum fluctu- 
ations extend over all length scales. Moreover, the amount 
of entanglement in quantum states is a valuable resource that 
promotes spin systems as candidates for quantum information 
devices 10,11 . 

Since the seminal studies on the interplay between en- 
tanglement and quantum critical fluctuations in spin 1/2 
models (j], l^t, several works suggested different local mea- 
sures of entangleme nt ( LME) as new tools to locate QCP's 
U H H EE El Hi- The term local here is meant for 
measures which depend on observables that are local in real 
space. This has to be contrasted withglobal measures of en- 
tanglement, e.g. the block entropy |2, 13], or the so-called 
localizable entanglement 1114111511 . which are aimed to capture 
the entanglement involved in many degrees of freedom. 

The picture that has emerged so far seems to be non sys- 
tematic and model-dependent. Some of the (local) indicators 
reach the maximum value at QCP's |9], while others show a 
singularity in their derivatives |5l l6l ll2ll . Close to the tran- 
sition, the system being more and more correlated, one ex- 
pects naively an increase of entanglement. However it seems 
that the maxima observed in single-site entropies have to be 
ascribed to a symmetry of the lattice Hamiltonian that does 
not necessarily correspond to a QCP. For example, in the ID 
Hubbard model the single-site entropy reaches the maximum 
possible value at U = [9]. This is due to the equipartition 
of the empty, singly-, and doubly-occupied sites rather than 
to the Berezinski-Kosterlitz-Thouless (BKT) transition occur- 
ring at that point. In fact, in the anisotropic spin-1 system 
discussed below, the equipartition points do not coincide with 
the transition lines which are not marked by any symmetry of 
the lattice model iflfjll . 

The onset of non-analyticity in two commonly used entan- 
glement indicators (concurrence and negativity) was recently 
proved in fnll for models with two-body interactions. Let us 



first argue, from a statistical mechanics point of view, that this 
result is in fact more general: as a consequence of the scal- 
ing hypothesis every local average displays a singularity at 
the transition with the exception of accidental cancellations. 
In particular as any LME is built upon a given reduced density 
matrix, the former will inherit the singularities of the entries 
of the latter. 

A second order quantum phase transition is characterized 
by long-ranged correlation functions and a diverging correla- 
tion length £. Let the transition be driven by a parameter g 
such that the Hamiltonian is 

U{g) = Ho + gV. 

At T — the free energy density reduces to the ground 
state energy density which shows a singularity in the second 
(or higher) derivatives with respect to g: 

-^(W(sO) = e(g) = e reg ( 5 ) + e sing (£(g)), 

where £ w \g— g c \ ~ v is the correlation length, g c is the critical 
point and L is the number of sites. Note that, as a consequence 
of the scaling hypothesis, the singular part of the energy e s i ng 
is a universal quantity that depends only on £, the relevant 
length scale close to the critical point. Hence, e S i ng may be 
considered quite in general an even function of g — g c around 
the critical point. 

Differentiating e(g) with respect to g, gives the mean value 
(V) / L, whose singular part O g behaves as 

O g *sgn(g-g c )\g-g c \ p . (1) 

Scaling and dimensional arguments imply that p = 
(d + C) v — 1 where d is the spatial dimensionality and £ is 
the dynamical exponent. For the sake of clarity here we set 
C = 1 as occurs in most cases fisil . For a second order phase 
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transition p > 0. In particular, if < p < 1 the next derivative 
will show a divergence 12511 
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In the case where g is mapped to the temperature T in the re- 
lated (d + l)-statistical model, C g will correspond to the spe- 
cific heat and p = 1 — a (Josephson's scaling law). As far 
as entanglement is concerned, the singular term O g appears 
in every reduced density matrix containing at least the sites 
connected by the operator V. Obviously, modulo accidental 
cancellations, any function (i.e. entanglement measures) de- 
pending on such density matrix, displays a singularity with 
an exponent related to p. The renormalization group theory 
allows us to be even more general: to the extent that a local 
operator can be expanded in terms of the scaling operators 
(permitted by the symmetries of the Hamiltonian), its average 
will show a singularity controlled by the scaling dimension of 
the most relevant term. 

From an operational point of view, LME's have been em- 
ployed mainly to detect the transition point using finite-size 
data. Following the previous discussion we can argue that, in a 
typical situation, the best suited operator for a finite-size scal- 
ing (FSS) analysis is precisely (V) for the following reasons. 
First because it naturally contains the most relevant operator, 
whose average O g has the smallest possible critical exponent 
p. Second the occurrence of sgn(g — g c ) in Eq. Q plays an 
important role in finding the critical point, in case its loca- 
tion is not known from analytical arguments. The FSS theory 
asserts 1 19] that in a system of length L, 



O g (L) ^ sgn(g - g c )L-p/^ o (L/0, 



(2) 



where $o( z ) is a universal function which must behave as 
z p l v in order to recover Eq. Q in the (off-critical) thermo- 
dynamic limit L ^> £. In the critical regime z — > 0, &o( z ) 
must vanish in order to avoid jump discontinuities for finite 
L. Notice that the sign of the microscopic driving parameter 
(.9 — 9c) survives in the FSS for O g . As a consequence, since 
O g (L) is an odd function of (g — g c ), the curves O g (L) at 
two successive values of L as a function of g cross at a single 
point g* L near g c (see below). In this way, by extrapolating the 
sequence g* L to L — > oo one has a useful method for detecting 
numerically the critical point. Surprisingly, to our knowledge 
the present method was not considered in the past, in favor of 
the so-called phenomenological renormalization group (PRG) 
method 1 19]. However the PRG method exploits the scaling 
of the finite-size gap which requires the additional calculation 
of an excited level typically computed with less accuracy than 
the ground state. This means that the computational time is 
roughly doubled. Another advantage w.r.t. the PRG method 
is that we do not have the complication of two crossing points 
g* L . In fact as A (g) ~ ~ \g — g c \ v is an even function, 
the scaled gaps will cross at two values of g for each L. 

Once g c is determined, using FSS techniques the critical 
exponents p and v may be extracted simply by estimating 



p/v = d + 1 — X/v, In order to find other possible critical 
exponents, we should perturb our model with other operators 
permitted by the symmetries H — > H + g'V' and repeat the 
same study near g' = 0. 

In what follows we will illustrate these ideas in two differ- 
ent d = 1 spin models: i) the thoroughly studied, exactly solv- 
able Ising model in transverse field for which all arguments 
can be checked analytically and ii) the spin-1 XX Z Heisen- 
berg chain with single-ion anisotropy for which there are no 
analytical methods to locate the different transition lines. 

Ising Model in Transverse Field. We consider the follow- 
ing Hamiltonian with periodic boundary conditions (PBC) 



H = - 
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where the a a, s are the Pauli matrices. This model exhibits 
a QCP at h = X, where it belongs to same universality class 
as the 2D classical Ising model, with central charge c = 1/2. 
From the exact solution, it is possible to show that the trans- 
verse magnetization m z — (erf), obtained differentiating the 
energy w.r.t. the driving parameter h, has the following ex- 
pression near the transition point h = X 
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As expected, m z is a continuous function at the transition 
point h = 1, showing a singular part which is manifestly odd 
in h — X. The next /i-derivative exhibits a logarithmic diver- 
gence, as it is related to the specific heat in the corresponding 
2D classical model. Most important for us is the "crossing 
effect" near the critical point of the family of curves m z (L), 
for different L. In Fig. [^to z (£) is plotted for several sys- 
tem sizes. It is evident that, increasing L, the crossing points 
converge rapidly to h = X. 

A quantitative analysis of the crossing effect may be done 
in the spirit of FSS, considering separately the critical regime 
(L <C and the off-critical one (L ^s> £) for any finite L. In 
the off-critical regime, the finiteness of the correlation length 
£, reflects in the exponential convergence of the energy to the 
thermodynamic value 
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with £ given by the formula sinh (l/2£) = |1 - h\ \h\ 1/2 /2, 
from which we can read the critical exponent v = 1, for h — > 
1 . In the critical regime, the finite-size expression for the 
transverse magnetization is 
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where -fc = 0.5772 ... is the Euler-Mascheroni constant. 
The finite-size critical field h c ,L is obtained via the cross- 
ing points between curves with slightly different lengths, 
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Figure 1: The transverse magnetization (erf), is plotted versus h for 
various sizes L ranging from 20 to 100 in steps of 10. The black thick 
line corresponds to the thermodynamic limit. The arrows indicate the 
direction of increasing L. The inset shows the derivative w.r.t. h. 



mL (h c ,L) = m t+2 (^c,l)- The solutionis 



h c ,L = 1 



showing a convergence towards the critical point as fast as 
L- 2 . 

As we stressed already, the singularities of local averages 
reflect in the behavior of LME's. Among these, the sim- 
plest measures the entanglement between one site and the 
rest of the system and is given by the von Neumann entropy 
Si = — Tipi In pi, where pi is the reduced single-site density 
matrix. For spin- 1/2 systems pi is simply written in terms of 
Pauli matrices 
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Pi 



(I + m x a x + m y (T y + m z a z ) 



For the Ising model m y — m x — and the single site 
entropy behaves as Si ~ —0.239 (h — 1) In \h — 1| so that 
its /i-derivative diverges logarithmically. A non-zero value 
of m x is possible if spontaneous symmetry breaking is taken 
into account, by adding a small longitudinal (i.e. along x) 
field that tends to zero after the thermodynamic limit is per- 
formed. In this case af becomes the most relevant operator 

1/8 

and m, = 9 (1 — h) (1 — /i 2 ) . Accordingly the singular 
part of the entropy is Si ~ (1 — h) 1 ^ 4 for h < 1. The same 
singularities are encountered in all the single-site measures 
built upon (pi) 2 , e.g. purity and linear entropy j6[. 

On the same line one can consider LME's based on the 
two-site density matrix pij, obtained taking the partial trace 
over all sites except i and j. The entries of pij now de- 
pend also on the two-point correlation functions (erf c^). In 
accordance with the general theory, all such averages be- 
have as (h — 1) In \h — 1| close to the critical point. In 



the case of nearest-neighbour sites this explains the log- 
arithmic divergence in the first derivative of the concur- 
rence C(l), as found in H. Instead the leading singu- 
larity in the next-nearest neighbour concurrence C (2) = 
[<of <+2> - ( CT M+2> + («+2> " 1] /2 turns out to be of 
the form (h — 1) In \h — 1| ||5|]. This is due to the accidental 
cancellation of the (h — 1) In \h — 1| terms contained in the 
correlators. 

Spin-1 Heisenberg Chain with Anisotropics. Let us now 
consider the non-integrable spin- 1 model 



H = J2 \SfSf + i 



qy qy 
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which shows a rich phase diagram 1 16). It is known that the 
transition line between the large-D phase (where the spins 
tend to lie in the icy-plane) and the Haldane phase (charac- 
terized by non zero string order parameters) is described by 
a conformal field theory with central charge c = 1 12011 . 
This means that the critical exponents change continuously 
along the critical line. For the detection of the c = 1 criti- 
cal line, the PRG I2II1 or the twisted-boundary method flrjl 
have been used in the literature. We have tested the finite-size 
crossing method outlined above, fixing A = 2.59 for which 
previous studies ensure p < 1 infill . The driving parame- 
ter being now D, the quantity to consider is de/dD which, 
by translational invariance reduces to ((5?)^) = O d . In 
Fig.[2]we plot Od versus D for various sizes L. The crossing 
points of the curves for subsequent values of L, determined 
by O d (D,L) — O d (D, L+ 10), converge rapidly to the crit- 
ical point D r = 2.294 consistently with the phase diagram 
reported in 1 16]. 

The effective theory in the continuum limit of the model @ 
around the c = 1 line reduces to the sine-Gordon Hamiltonian 
density Eoll 
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The coefficient p is zero along the critical line, a is a short 
distance cut-off of the order of the lattice spacing and K is 
related to the compactification radius, varying continuously 
between 1/2 and 2 along the critical line. In this framework, 
crossing the critical line in the lattice model means going 
from negative to positive values of p and the corresponding 

/U-derivative gives C M = ^cos (^/AnK^j ^. From the sine- 
Gordon theory 1 23] it is known that C M « sgn (p) \p\ K I( 2 ~ K ^ 
and £ ~ \p^^ K ~ 2 \ In our case p w (D — D c ) at fixed 
A, so that p = K/(2-K),u = l/{2-K). On the one 
hand, the critical exponent pjv — K, can be independently 
calculated from the conformal spectrum obtained numerically, 
as explained in Ref. 1 20] giving K = 0.76. On the other hand, 
from the FSS of the derivatives of Od at D = D c (shown in 
the inset of Fig. [2) we find K — 0.78 showing that the method 
presented here is effective for the calculation of the critical 



4 



0.65 




Figure 2: The single-site average Od = {(S?) 2 ), is plotted versus D 
for various sizes L ranging from 20 to 100 (thick line) in steps of 10. 
The arrows indicate the direction of increasing L and the inset shows 
numerical derivative w.r.t. D, interpolated with splines. The data 
have been obtained via a DMRG program 12211 using 400 optimized 
states and 3 finite system iterations with PBC 



such singularities can be traced back to the behavior of the 
transition-driving term V and to the corresponding scaling di- 
mension. Moreover the FSS of (V) turns out to be a valuable 
method to determine the critical point and the associated ex- 
ponents. This method has been illustrated for a couple of spin 
models displaying qualitatively different QCP's. More gener- 
ally, these considerations can be directly transposed to other 
many-body problems, like strongly interacting fermionic sys- 
tems. Our arguments indicate that the singular behavior of 
LME's can be adequately understood in terms of statistical- 
mechanics concepts. Physically, the understanding of the inti- 
mate relation between genuine multipartite entanglement and 
the critical state remains an open challenge. From this per- 
spective, it may be useful to conceive nonlocal indicators that 
could unveil the role of non-classical correlations near criti- 
cality. 

We are grateful to E. Ercolessi, F. Ortolani and S. Pasini for 
useful discussions. This work was supported by the TMR net- 
work EUCLID (No. HPRN-CT-2002-00325), and the COFIN 
projects 2002024522.001 and 2003029498.013. 



exponent as well. Since the same transition can be driven by A 
at fixed D, we checked that sitting at D = 2.294 we obtained 
A c = 2.591 by looking at Ox = (SfSf +1 ). According to out- 
general discussion, a similar behavior is seen also in Od even 
if it is a single-site indicator. 

The scaling exponent p/v in Eq. is best obtained from 
the analysis of the first derivative when p < 1. As we move 
towards the BKT point, K — > 2, v — > ac, so the divergence 
should be seeked in derivatives with increasing order. Ac- 
cordingly the crossing method (as well as the PRG) becomes 
less efficient as we approach the BKT point, for which a finer 
analysis is needed involving level spectroscopy p24j- 

Again the singularities of local averages enter the LME's. 
In the spin-1 case, thanks to the symmetries of the Hamilto- 
nian Q the the single-site entropy reads 

5i = -O d log f^j - (1 - O d ) log (1 - O d ) , 

where < Od < 1 in any bounded region of the phase di- 
agram. Note that the maximum of Si occurs for Od = 2 /3, 
which is not related to any phase transition, but simply signals 
the equipartition between the three states | + 1), |0), | — 1). 
This occurs for example at the isotropic point (A = 1, D = 0) 
where the system is known to be gapped. Similarly, in the 
Ising model Si is maximal at h = 0, i.e. when m z = 0, 
where no transition occurs. Therefore the intuitive idea of the 
local entropy Si being maximal as a criterion to find quantum 
phase transitions |9], seems to be more related to symmetry 
arguments rather than to criticality. 

In this Letter we have put in evidence the origin of singu- 
larities in LME's which have been recently proposed to de- 
tect QCP's. Typically, apart from accidental cancellations, 
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